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We report the results of some recent E791 charm analyses. They include 1) a search for rare 
and forbidden decays, 2) measurements of form factors for D + —* K*°£ + u e and D$ — > (f>£ + v e , and 
3) Dg and D° lifetime measurements including the lifetime difference between D° — > K~ir + and 
D° — » K~ K + . The latter is the first direct search for a possible lifetime difference that could 
contribute to D° — D mixing. 
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I. INTRODUCTION 

as : 
as . 

Q\ • E791 is a high statistics charm experiment that acquired data at Fermilab during the 1991-1992 fixed-target run. 
The experiment combined a fast data acquisition system with an open trigger. Over 2 x 10 10 events were collected 
with the Tagged Photon Spectrometer |jj using a 500 GeV tt~ beam. There were five target foils with 15 mm center- 
<^ to-center separations: one 0.5 mm thick platinum foil followed by four 1.6 mm thick diamond foils. The spectrometer 
■ included 23 planes of silicon microstrip detectors (6 upstream and 17 downstream of the target), 2 dipole magnets, 
I , 10 planes of proportional wire chambers (8 upstream and 2 downstream of the target), 35 drift chamber planes, 2 
multi-cell Cerenkov counters that provided ir/K separation in the 6-60 GeV/c momentum range Q, electromagnetic 
' and hadronic calorimeters, and a muon detector. 

i> : 

SO ■ 

O ' II. RARE AND FORBIDDEN DECAYS OF D+ , D s AND D° 

o ■ 

Q\ ■ We report the preliminary results of a search for flavor-changing neutral-current (FCNC), lepton number violating 
(LNV) and lepton family violating (LFV) decays of D + , D$ and D°, into modes[| containing muons and electrons. 
^ | This analysis is an extension of our previous work [|| . The search for rare and forbidden decays provides another test 
^ ■ of the Standard Model. The contribution, within the Standard Model, of short distance electroweak processes to the 
^H. branching ratios of the FCNC-like decay modes described here is expected Q to be less than 10 -8 . The LNV and 
LFV decay modes are totally forbidden by the Standard Model. Therefore, the search for these rare and forbidden 
decay modes allows tests for violations of the Standard Model from FCNC, neutrino oscillations or other even more 
' exotic processes. We are currently examining all the modes described in Table ^ however, we are only presenting the 
following searches in this paper: D + —* n + e + e~ , D + — > 7r + /i ± e =F , D + — > 7r~e + e + , D + — > 7r~fi + e + , D° — > e + e~, and 

TABLE I. Rare and Forbidden Decay Modes. 



FCNC 


LFV 


LNV 


D + -f iv+n+fj,- 


D+ -f TT+^e^ 


D + -> 7T-/X+/X+ 


D+ ~* iv+e+e- 


D+ -» 7r-/i + e + 


D+ -> Tv-e+e + 


D+ K+^+fi- 


D+ -> K+t^e^ 


D+ K-^+^+ 


Dg -> K+e+e- 


D+ -> K-^+e+ 


Dt -> K~e + e + 


D° -» p+p- 


D° -► ^e* 





D° 



1 all references to D + , D$ and D° and their decay modes also imply the corresponding charge-conjugate states. 
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To separate charm candidates from background, we require the secondary vertex to be well-separated from the 
primary vertex and located well outside the target foils and other solid material, the momentum vector of the candidate 
D + or D° to point back to the primary vertex, and the decay track candidates to pass closer to the secondary vertex 
than the primary vertex. Specifically, the secondary vertex must be separated by > 20 a l for D + , (> 12 a l for D°) 
from the primary vertex and by > 5.0 from the closest material in the target foils, where <7j, in each case is the 
calculated longitudinal resolution in the measured separation. The sum of the momentum vectors of the three tracks 
from this secondary vertex must miss the primary vertex by < 40 /xm in the plane perpendicular to the beam. We 
form the ratio of each track's smallest distance from the secondary vertex to its smallest distance from the primary 
vertex, and require the product of these ratios to be less than 10~™ where n is the number of tracks in the secondary 
vertex. There is also a cut on the lifetime of the D + or D° candidate; it is < 5 picoseconds for the D + and < 3 
picoseconds for the D . Finally, the net momentum of the D + candidate transverse to the line connecting the primary 
and secondary vertices must be less than 250 MeV/c (300 MeV/c for the D°). 

For this analysis we use a "blind" analysis technique. This method entailed the following steps: first covering the 
signal region with a "box" , then optimizing ALL cuts using Monte Carlo simulated signal events and data from the 
wings outside the "box" before opening the "box" and only then opening the "box" that covers the signal region. The 
closed "box" mass windows are: 1.84 < M{D+) < 1.90 GeV/c 2 for D+ -> nfjtfi, 1.78 < M(D+) < 1.90 GeV/c 2 for 
D+ -> Tree and nfie, 1.95 < M(D$) < 1.99 GeV/c 2 for -> Kfifj,, 1.91 < M(D$) < 1.99 GeV/c 2 for -> Kee 
and Kfie, 1.83 < M(D°) < 1.90 GeV/c 2 for D° -> fifi, and 1.76 < M(D°) < 1.90 GeV/c 2 for D° -> ee and fie. One 
should note the asymmetric windows for the decay modes containing electrons. This is to account for the electron 
bremsstrahlung tail. The upper limit of the branching ratio is calculated using the following formulae: 

RR Nx/sx Rp _ N x e norm 

OKx — -TZ j • Dtinorm — "TT ' ^ ti n orm Uj 

'norm/ ^-norva ^ norm £-X 

where Nx is the number of observed events for a given decay mode and e is the detection efficiency. The ratio of the 
detector efficiencies for the normalization decay mode versus the unknown decay mode is 

N MC 



ex ' 7VMC 



Nx is corrected, using the method of Feldman and Cousins ||, for background. We have not yet corrected for 
systematic errors. The number of normalization events for the D + decay modes is N norm = 25330 ± 169, D + — > 
K~ir + ir + events, and for the D° decay modes is N norm — 26980 ± 185, D a — > K~ir + events. 

In addition to the cuts made on the kinematic variables described above, particle identification cuts are used to 
reduce the background. These cuts tagged the dilepton candidates, thereby greatly reducing the background. Events 
which could also be reconstructed as a reflection, due to particle misidentification, from other purely hadronic decay 
modes into a region outside the "boxes" , but within the overall mass window, were removed. Examples of this would 
include D + — > K~tt + tt + reconstructed as D + — > n~£ + £ + or jjo _> K~tt + reconstructed as D° — > £ + £~. In Figure [|, 
we present preliminary plots of the data for the D + and D decay modes described in this paper. The yellow boxes 
are the part of the plot that were originally covered by the "boxes" . The events observed within these "boxes" are 
the total number of observed events, signal plus background. In the plots of the D + decay modes the lower yellow 
box is the D + "box" and the upper yellow box is the Dg "box" . 
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FIG. 1. The preliminary data for the D + and D° decay modes presented here. 
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We estimate the mean background assuming a flat distribution calculated from the events seen in the regions that 
are outside the "boxes" . Then the Poisson signal mean fj, for the total number of observed events for a 90% confidence 
level upper limit are calculated, using the method of Feldman and Cousins [^), and this number, Nx, is substituted 
into Equation [j] thus giving a 90% confidence level upper limit branching ratio. The estimated background, total 
number of events within the "box" , the number N\ for a 90% confidence level upper limit, the calculated upper limit 
branching ratios and the previously measured upper limits Q are given in Table |. 
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TABLE II. Preliminary Results 90% C.L. Upper Limits 



Mode 


Estimated Background 


Observed 


Calculated Nx 


Branching Ratio 


1998 PDG BR 


D+ -> vr+e+e" 


0.6 


1 


3.76 


< 4.3 x 10 -5 


< 6.6 x 10" 5 


D + _> Tr+^igT 


0.0 


1 


4.36 


< 3.2 x 10~ 5 


< 1.2 x 10 -4 


D + -> 7r~e+e+ 


0.0 


2 


5.91 


< 7.7 x 10~ 5 


< 1.1 x 10~ 4 


D+ -» 7r-/i+e+ 


0.6 


1 


3.76 


< 3.6 x 10" 5 


< 1.1 x 10" 4 


D° -> e+e" 


1.8 





1.30 


< 5.2 x 10" 6 


< 1.3 x 10" 5 


£>° -> /i±e+ 


2.6 


3 


4.80 


< 1.0 x 10~ 5 


< 1.9 x 10~ 5 



a By Method of Feldman and Cousins 



A comparison between our calculated 90% confidence level upper limit branching ratios and the previous results 
is shown by the plot of the 90% C.L. upper limits in Figure g. 
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FIG. 2. The 90% C.L. upper limit branching ratios for E791 data and for the 1998 PDG values. 
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III. D s LIFETIME MEASUREMENTS 



Accurate measurements of the lifetimes of the weakly decaying charm mesons are useful for understanding the 
contributions of various weak decay mechanisms. Despite the fact that they are all tied to the charm quark decay, 
the decays of the ground state charm mesons can have different contributions from the four first-order processes (two 
spectator, W-annihilation, and W-exchange) and their lifetimes are, in fact, quite varied l|. 

r(P+) :r(P°) :r(P s ) = 2.5: 1:1.1 (3) 

We have made a new precise measurement of the Ds lifetime using 1662 ± 56 fully reconstructed Dg — > (fm + 
decays. The Ds lifetime is measured using an unbinned maximum-likelihood fit to be 0.518 ± 0.014 ± 0.007 ps. This 
value is somewhat higher than the world average D° lifetime j|] of 0.467 ± 0.017 ps. Using our result and the world 
average D° lifetime (§, we find that the ratio of the Ds lifetime to the D° lifetime to be 

r - = 1.25 ± 0.04 (a 6cr difference from unity) (4) 

showing, for the first time, significantly different lifetimes for the Ds and D°. This result may be used to constrain 
the contributions of various decay mechanisms to charm decay and further refine our quantitative understanding of 
the hierarchy of charm particle lifetimes. 



IV. MEASUREMENT OF THE LIFETIME DIFFERENCE BETWEEN D° -> K~ir + AND D° -» K~K+ 



We report the first directly measured constraint on the decay-width difference Ar for the mass eigenstates of the 
D° — D system ||. We obtain our result from lifetime measurements of the decays D° — » K~n + and D° — > K~K + , 
under the assumption of CP invariance, which implies that the CP eigenstates and the mass eigenstates are the 
same. The lifetime of the CP-even final state D° — > K~K + , as calculated from 6683 ± 161 weighted events, is 
tkk = 0.410 ± 0.011 ± 0.006 ps. The lifetime of an equal mixture of CP-odd and CP-even final states D° — » K~n + , 
as calculated from 60581 ±353 weighted events, is t K7I = 0.413 ± 0.003 ± 0.004 ps. We find that AT = 2 (T KK - T K7T ) 
= 0.04 ± 0.14 ± 0.05 ps" 1 , leading to a limit of -0.20 < Ar < 0.28 ps" 1 at 90% C.L. The value of Ar is consistent 
with zero and thus, at our level of sensitivity, is consistent with the Standard Model. 



V. MEASUREMENTS OF D + -> K*°£+u t AND D+ -> <j>t+v t FORM FACTORS 



The Vector and Axial form factors are used to extract elements of the Cabibbo-Kobayashi-Maskawa mixing matrix. 
Here we present measurements of the form factor ratios, evaluated at q 2 = (GeV/c) 2 , ry — V (0) j A\ (0), r-i = 
A2 (0) fA\ (0), for both the muon and electron channels, and r3 = A3 (0) /A% (0) in the muon channel. This is the 
first simultaneous measurement of both the muon and electron channels. 



A. D + -> K*°e+u e Form Factor Ratios 



We have already reported the form factor ratios ry and r2 for the decay D + — > K*°e + v e , K*° — > K 7r + [g| . The 
form factor ratios in the decay D + — > K*° ^, + v^, K*° —> K~ir + were measured using 3034(595) signal (background) 
events and yield §r F = 1.84±0.11 ±0.09, r 2 = 0.75±0.08±0.09, and as a first measurement r 3 = 0.04±0.33±0.29. 
The combined electron and muon form factor ratios, for approximately 6000 events, are ry = 1.87 ± 0.08 ± 0.07 and 
r 2 =0.73 ±0.06 ±0.08. 
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B. — ► (f>£ + P e Form Factor Ratios 



We have measured the form factor ratios ry and for the decay Dg — > 4>l + v^ tfi — ► K~ K + in both the electron and 
muon channels |Tjj| . The form factor ratios in the electron channel were measured using 144(22) signal (background) 
events and yield ry = 2.24 ± 0.47 ± 0.21 and r-i = 1.64 ± 0.34 ± 0.20. In the muon channel there were 127(34) signal 
(background) events, yielding ry — 2.32 ± 0.54 ± 0.26 and r2 = 1.49 ± 0.36 ± 0.20. The combined electron and muon 
form factor ratios, for approximately 271 events, are ry = 2.27 ± 0.35 ± 0.22 and — 1.57 ± 0.25 ± 0.19. 
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